Advantage of U+U over Au+Au collisions at constant beam energy 
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Collisions of deformed uranium nuclei are studied in a Monte-Carlo Glauber model. For U+U at 
zero impact parameter (b = 0) in the most favorable orientation (tip-to-tip), the transverse particle 
density (charged-particle rapidity density per weighted transverse area of the initial participant 
zone) increases by about 35% compared to Au+Au at b = 0. To estimate the advantage of U+U 
over Au+Au in the context of real experiments at the Relativistic Heavy Ion Collider, we examine 
the effect of a range of centrality cuts on the event sample. In terms of the transverse particle 
density, the predicted advantage of U+U is about 16%. 

PACS numbers: 25.75.-q,25.75.Ld,24.85.+p 



I. INTRODUCTION 



Large space-time volumes with the highest energy den- 
sity to date, over 15 GeV/fm 3 in typical estimates 0, 
have been produced in ^snn = 200 GeV Au+Au col- 
lisions at the Relativistic Heavy Ion Collider (RHIC) 
Hll 110 At such high energy densities, strongly in- 
teracting matter is predicted to be in a phase commonly 
referred to as quark-gluon plasma (QGP) A goal of 
the RHIC program is to study this phase in the extended 
systems created in relativistic collisions of heavy nuclei. 

The presence of QGP manifests itself in the equation of 
state, which (together with the initial conditions) deter- 
mines various measured properties and serves as an input 
to fluid-dynamical calculations. It is reasonable to expect 
that, barring calculational problems, the main features of 
the equation of state can be probed by comparing fluid- 
dynamically calculated quantities to data. It has been 
argued recently that the quark-gluon matter produced 
in 200 GeV central Au+Au collisions behaves as a per- 
fect fluid with negligible viscosity This conclusion 
was based in part on the agreement [5| of the flow quan- 
tity vift (see Eqs. an d @ f° r a definition) in central 
Au+Au with the prediction of an ideal-fluid calculation 
@ . The data seem to indicate that the systems produced 
in 200 GeV central Au+Au collisions have just reached 
a large enough energy density in a sufficient space-time 
volume for the approximate agreement with the fluid- 
dynamical results. If energy density could be increased 
even further, a crucial test would be whether «2/e con- 
tinues to grow, or whether it saturates at the value for 
an ideal fluid. 

A promising avenue to increase the energy density 
and/or increase the volume of high-density matter with- 
out any increase in beam energy is to collide heavy de- 
formed (prolate) nuclei. This has been proposed by sev- 
eral theorists [ifi fill Il2l| . Central collisions with the long 
axes aligned with the beam (we call this the tip-to-tip 
configuration, as in Ref. ^lj) would be the most desired 
configuration. However, as long as no beams of aligned 
deformed nuclei are available [lJl, the desired configura- 
tions have to be selected bv exDerimenters usine: some 



combination of triggering and offline event selection. Re- 
cently, Heinz and Kuhlman advocated the use of full- 
overlap collisions between deformed uranium nuclei |l2j . 
The selection of full-overlap collisions using Zero-Degree 
Calorimeters (ZDCs) was examined in Ref. 01- How- 
ever, various fluctuations and background signals in the 
detectors mean that an event sample with the lowest ZDC 
signal will still include a substantial fraction of collisions 
in which full overlap did not occur. In this work, we ad- 
dress these issues with the help of data on the relevant 
performance of the STAR detector [l5j . 

The paper is organized as follows: in Scc.[H]wc briefly 
review the relevant quantities and outline our calculation, 
while in Sec. IIIII we present results without and with 
detector smearing. In Sec. IIVI we briefly summarize our 
findings. 



II. BACKGROUND AND CALCULATIONAL 
FRAMEWORK 

The elliptic flow anisotropy 1)2, defined as the second 
Fourier coefficient in the expansion 
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of the azimuthal distribution dN/d& of final state par- 
ticles is a sensitive measure of the success of fluid- 
dynamical models. Based on suggestions in Refs. 0, 
IT?! ITsj it has become customary to consider v 2 /e, the 
elliptic flow scaled by the initial spatial eccentricity 



Or 



[u 2 



(2) 



The normalization by e emphasizes that the final momen- 
tum anisotropy is driven by the initial geometry. One 
advantage of using V2 /e is that the impact-parameter de- 
pendence is largely removed from its fluid-dynamically 
calculated value. This quantity is frequently plotted 
against the transverse particle density (1/S) dN c h/dy, 
where dN„h I dv is the multiolicitv of chareed Darticles 
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per unit rapidity and 



s = tt V(^ 2 > (y 2 ) 
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is the transverse area of the overla p zo ne weighted by 
the number of wounded nucleons 0, El • The trans- 
verse particle density (1/5) dN^/dy can be interpreted 
as a measure of the initial entropy density in the trans- 
verse plane [l2T |. An increase in this quantity means in- 
creased particle production, and is thus associated with 
an increase in energy density in the system. The charged 
particle multiplicity per unit rapidity dN c h/dy can be 
increased by increasing the number of binary collisions, 
achievable by tip-to-tip collisions of heavy deformed nu- 
clei. Furthermore, the U+U tip-to-tip configuration has 
about the same overlap area S as Au+Au collisions at 
6 = 0, and thus contributes to a further increase in 
(l/iS) dN^/dy beyond what would be obtained with hy- 
pothetical spherical nuclei of the same mass as uranium. 
On the other hand, for b = collisions, S increases by 
~ 24% between Au+Au and U+U when the latter are 
in the body-to-body orientation, i.e., when the long axes 
of the nuclei are parallel to each other, but perpendic- 
ular to the beam. Thus it can be concluded that even 
in the experimentally unrealistic limit where full-overlap 
U+U collisions can be isolated, it is necessary to also 
distinguish between different full-overlap configurations 
in order to realize the full potential advantage of U+U 
collisions. 

The RHIC data as a function of (1/-S) dN c h/dy reach 
the perfect-fluid value of i^/e ~ 0.2 in the most central 
200 GeV Au+Au collisions at around (1/ 'S)dN ch /dy ~ 
28 fm~ 2 . Since viscosity would decrease the predicted 
Vi/e, it would be difficult for a fluid-dynamical model to 
be reconciled with any further increase in the observed 
ratio Vijt. Thus, the behavior of this ratio at signifi- 
cantly larger (1/ S) dN c h/dy should prove informative in 
assessing the applicability of fluid dynamics in central 
collisions at RHIC. 

Uranium ( 238 U) beams have been proposed as a way 
to increase (1/ 'S)dN c h/dy 0, 0- We assume in 
this study that, in order to isolate a sample of U+U col- 
lisions with the maximum possible (l/S) dN c h/dy, the 
best available experimental procedure is to select events 
where the fewest possible spectator neutrons are detected 
(low ZDC signal) in coincidence with the highest ob- 
served multiplicity of produced particles in a broad re- 
gion centered on mid-rapidity. Available data from the 
STAR experiment ^5| describe this correlation and its 
fluctuation from event to event, and allow the output 
from the Monte-Carlo Glauber model described below to 
be filtered appropriately to simulate the relevant experi- 
mental limitations. 

We represent the quadrupole deformation of the 
ground-state uranium nucleus in the standard |20| way: 
we take a Saxon- Woods density distribution with sur- 
face thickness a = 0.535 fm and with R = i? sp (0.91 + 
0.27 cos 2 9), where 9 is the polar angle relative to the sym- 
metry axis of the nucleus and R SD = 1.12A 1 / 3 -0.86A" 1 / 3 



fm |20j . The small hexadecapole moment of the uranium 
nucleus is neglected. This yields i?iong/-R P crp — 1-29 
|9l lTolll2T |. The orientation of the first and second nucleus 
in the colliding pair is fixed by the two angles (9 p , <f) p ) 
and (9 t , 0t), respectively. The angles 9 p and 9 t describe 
the orientation of the symmetry axis relative to the beam 
direction, and they are uniformly distributed in [0, 7r/2]. 
The azimuthal angles </> p and 4>t describe rotations about 
the beam direction, and they are uniformly distributed in 
[0, 27r]. A schematic view of a cut in the transverse plane 
illustrates the collision geometry in Fig.^ where the lines 
represent equivalent sharp surfaces, and the shaded area 
corresponds to the overlap region in a collision at im- 
pact parameter b with fixed orientations of both nuclei. 
To simulate a uranium nucleus, 238 nucleons are selected 
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FIG. 1: (Color online) Collision geometry. 



randomly according to the appropriate Saxon- Woods dis- 
tributions with distances between any two nucleons satis- 
fying r.ij > 0.4 fm. In the collision process, nucleons arc 
considered wounded if the transverse distance between 
them becomes < y/ A.2/ir fm, where we use 4.2 fm 2 for 
the nucleon-nucleon cross section at 200 GeV. To convert 
between track densities in rapidity y and in pseudorapid- 
ity 77, we use the approximation dN c h/dy ~ 1.15 dN c h/dr] 
and apply the parameterization 



dNch 
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Tlpp [xN b + (1 - x)N w /2] 



(4) 



where Nb is the numbers of binary collisions and N w is 
the number of wounded nucleons \22\ . The values of the 
parameters used are x = 0.15 and n pp = 2.19 at 200 GeV. 
These values provide a reasonable fit to the PHOBOS 
data |2^| as shown in Fig. [2] 



III. RESULTS WITHOUT AND WITH 
DETECTOR SMEARING 

A. No smearing 

Figure |3 displays the charged multiplicity distribution 
dN c h/dri for U+U collisions. Since we average over all 
possible orientations, this distribution is similar to that 
obtained for hypothetical spherical nuclei with mass 238. 
For ideal tip-to-tip U+U collisions (6 = fm, 9 D = 0, 9 t 
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PHOBOS data Au + Au at 200 GeV 
MC Glauber Au + Au at 200 GeV 
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FIG. 2: (Color online) Fit to PHOBOS data for Au+Au at 
200 GeV, using Eq. gj. 




400 600 800 1000 1200 
dN ch /dr, 

FIG. 3: Charged particle multiplicity distribution for U+U 
collisions. 



= 0) wc obtain (l/S)dN ch /dy ~ 42.6 fm~ 2 . For body- 
to-body collisions (b = 0, 9 P — 9t = 7r/2, <t> P = 4>t = 0) 
we have (l/S)dN ch /dy ~ 31.7 fm -2 , surprisingly close 
to 31.5 fm -2 , the value for Au+Au at b = fm. 

Although dN c h/dy is larger for U+U than for Au+Au, 
the larger overlap area for body-to-body collisions offsets 
the increase and results in a value of (1/S) dN c h/dy close 
to that for Au+Au. The approximately 35% increase 
in (1/S) dN c h/dy in tip-to-tip U+U collisions compared 
to Au+Au is due to the increase in dN c h/dy, since the 
overlap areas are similar in these cases. 

Under the best of circumstances, an experiment can 
hope to select certain ranges of angles and impact pa- 
rameter. For illustrative purposes, we choose a somewhat 
arbitrary range of definitions for near tip-to-tip configu- 
rations: b < 1 or 2 fm, with 9 p and 9 t both within 10° or 
20°, and for simplicity we keep 8 P = 9t throughout. 

Even if we impose the strictest definition above for 
near tip-to-tip events (9 p = 9 t < 10° and b < 1 fm) 
and a similarly restrictive definition for near body-to- 
body collisions, almost all of the near tip-to-tip and near 
body-to-body events lie in the top 3% of the dN c h/dr] 
(or, equivalcntly dNch/dy) distribution. For Au+Au col- 



lisions with dN c h/dy in the top 3% of the event sample, 
the additional requirement that the number of spectator 
nucleons lie in the bottom 1% of its parent distribution 
leads to little or no further change in (1/S) dN c h/dy. For 
near tip-to-tip U+U configurations, there are significant 
increases of 32-36% over central Au+Au. Details related 
to various near tip-to-tip U+U selection conditions can 
be found in Table I, and as before, a spectator cut has 
negligible additional effect on (1/ S) dN c h/dy. If a near 
body-to-body configuration is defined as < b < 1 fm, 
9 P , 9t > 80°, and \<j> p — 4>t\ < 10°, then the mean value of 
(1/5) dN c h/dy comes out to be 31.3 fm -2 , very close to 
the ideal body-to-body case mentioned above, and still 
negligibly different from the value for central Au+Au. 





spectator 


< b < 1 fm 


< b < 2 fm 




count 


%evt 


/ 1 dN ah \ 
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%evt 


/ 1 dN ah \ 
\S dy 1 


e < io° 


any 


0.174 


42.25 


0.687 


41.59 


6 < 20° 


any 


0.707 


41.75 


2.756 


41.13 


e < io° 


lowest 1% 


0.173 


42.26 


0.478 


41.95 


6 < 20° 


lowest 1% 


0.704 


41.76 


1.903 


41.50 



TABLE I: Mean values of (1/S) dN c h/dy and percentage of 
the event class in near tip-to-tip U+U configurations relative 
to the number of events in the top 3% of dN c h/dy. 

In the discussion above, central collisions were selected 
primarily via their high multiplicity. Next, we use a low 
spectator count as the initial centrality criterion, and ex- 
plore the effect of an additional cut on total multiplicity. 
Results are presented in Fig. 0] and it is evident that the 
ratio of near tip-to-tip to near body-to-body events in- 
creases very substantially as more central collisions are 
selected. However, even for the 5% of collisions with 
highest multiplicity, the near tip-to-tip category remains 
only about 21% of the events with a low spectator count 
(the latter being defined as the lowest 1% of all collisions, 
regardless of multiplicity). 



B. With smearing 

In the previous section, we defined ranges of U+U col- 
lision configurations termed "near tip-to-tip" and "near 
body-to-body" . These ranges were defined for illustrative 
purposes only; in this section, we estimate the character- 
istics of a class of events intended to have the highest 
experimentally obtainable enrichment of the desired tip- 
to-tip configurations. We assume that this is carried out 
by selecting events with the highest charged particle mul- 
tiplicity in a broad rj region centered on midrapidity, and 
simultaneously cutting to select events with a low signal 
in both Zero Degree Calorimeters. 

The measured dN c h/drj in a detector like STAR is lin- 
early correlated with the true dN c h/drj, but is subject 
to random event-to-event fluctuations in the efficiency 
for detecting tracks as well as random variations in the 
inclusion rate for tracks that do not originate from the 
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FIG. 4: Among events where the spectator micleon count 
lies in the lowest 1% of all collisions, the percentage of near 
tip-to-tip (6 < 20°, b < 2 fm) and near body-to-body 
(6 > 70°, \<j)p - <t>t\ < 20°, & < 2 fm) configurations as a 
function of an additional cut on dN c h/dr]. The horizontal 
scale is normalized such that 100% corresponds to the full 
sample where only the spectator cut has been applied. 



primary collision vertex. Similarly, the ZDC signals have 
a component that is proportional to the number of emit- 
ted spectator neutrons as well as a component from var- 
ious sources of background. Furthermore, the RHIC 
intersection regions arc arranged such that the ZDCs 
have acceptance for only free neutrons and an insignif- 
icantly small fraction of neutron-rich charged spectator 
fragments. Thus while we compute the total number of 
spectator nucleons in a given event using our Monte Carlo 
Glauber model, in real collisions at STAR, only a subset 
of the spectators are detected by the ZDCs, and so events 
with a small ZDC signal include a high background level 
of collisions with incomplete overlap. 

The upper panel of Fig.[S]shows contours of event den- 
sity in the plane of spectator count versus dNch/dr] from 
our Monte Carlo Glauber model without any detector 
simulation. The resulting correlation is much narrower 
than what is observed in experiment [l5j . an expected 
outcome given the simplicity of the model and the neglect 
of the experimental effects discussed above. To make 
a realistic estimate of the best enrichment of tip-to-tip 
U+U collisions that might be achieved, it is necessary to 
smear out the narrow ridge in the upper panel of Fig. [S] 
so that a pattern resembling the experimental distribu- 
tion is obtained. For this purpose, a Gaussian-distributed 
random number is first added to dN c h/drj for each event, 
such that the steeply dropping upper tail of the dN c h/dr] 
distribution in Fig. |3| is smeared out to the point where 
it resembles the experimental distribution from STAR. 
The smeared multiplicity distribution is shown in Fig. [5] 
Note that we are concerned here with the upper end of 
the multiplicity distribution (highest values of dN c h/dr]), 
and the variation of the smearing at lower multiplicities 
is not of interest. 

Next, another Gaussian-distributed random number 
(with non-zero mean) is added to the spectator count 



and the result is used as the simulated ZDC signal. The 
Gaussian mean and width are adjusted to produce a ridge 
near the lower right part of the lower panel of Fig. [5] that 
resembles the observed data from STAR . This proce- 
dure has a small effect on the agreement of the calculation 
with the PHOBOS data shown in Fig. [21 and while the 
effect is minimal, we corrected for the difference. 
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FIG. 5: Contours of event density in the plane of simu- 
lated ZDC signal versus dN c h/drj before (upper panel) and 
after (lower panel) simulation of detector-related smearing 
and background. The contour levels are the same in both 
panels, with a ratio of 1:7 between the two outermost con- 
tours, and a ratio of 1:2 between all other adjacent contours. 



Various combinations of cuts in dN c h/ drj and in ZDC 
signal have been applied to the simulated U+U data plot- 
ted in the lower panel of Fig.[5]with the objective of max- 
imizing the transverse particle density (1/ S) dN c h/dy. 
Because the quantity we seek to maximize is trivially 
correlated with one of our cut variables, we assess the 
maximum particle density in U+U collisions relative to 
Au+Au collisions with the same cuts on a percentage 
basis. 

The vertical axis in Fig. [7| shows the ratio 
[(l/S) dN ch /dy] V u/[(l/S) dN ch /dy] AuAu as a function of 
independent scans across both cut variables. Note the 
zero-suppressed vertical axis; we conclude that no fur- 
ther increase in [(1/ S) dN c h/dy]\jjj beyond about 18% 
relative to Au+Au can be reached in the context of our 
simulation, even searching beyond the plotted combina- 
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continuing to the highest RHIC energy. On the same 
number line, we indicate that when our simulation for 
Au+Au is constrained to agree well with the value from 
the STAR collaboration, we predict a potential further 
increase of about 4 fm -2 for U+U with the same central- 
ity cut. The calculated value of (1/5) dN c h/dy is larger 
than the STAR result for Au+Au at 200 GeV because of 
different parameterization of the nuclear radius and small 
uncertainties in the number of wounded nucleons. How- 
ever, this has a negligible effect on the relative increase 
from Au+Au to U+U. 



FIG. 6: Charged particle multiplicity distribution in U+U 
collisions after detector smearing. 




FIG. 7: Ratio [(1/5) dN ch /dy} vv /[(1/S) dN ch /dy] A uAu as a 
function of different dN c h/dr) and ZDC cuts, labeled according 
to the selected percentage of the total sample. 



tions of cuts on dN c h/dr] and ZDC signal. In the region 
of the most selective combination of centrality cuts, the 
advantage of U+U is about 16%. 
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FIG. 8: Maximum values of (1 / S) dN c h / dy . Values from 
NA49 and STAR correspond to central Pb+Pb [l<J and cen- 
tral Au+Au H, respectively. The results labeled Au+Au 
ysTw = 200 GeV and U+U v /JKw = 200 GeV correspond 
to the top 5% in dN c h/drj from the present simulation. 



Figure [S] summarizes the measured progression in the 
maximum {l/S)dN c h/dy for central Pb+Pb or Au+Au 
collisions, beginning with a low SPS beam energy and 



IV. CONCLUSION 

We have undertaken a study of the potential for U+U 
collisions to realize an increased transverse particle den- 
sity (\ / S)dN c h / dy without increasing the beam energy. 
We note that the present investigation addresses only a 
limited aspect of U+U collisions, namely, we have stud- 
ied the "worst case scenario" where only the transverse 
particle density (l/S)dN c h/dy is of interest, and where 
the charged particle multiplicity and the Zero-Degree 
Calorimeter signals are the only means of selecting the 
desired U+U events. Previous studies have considered 
the issue of the path-dependent energy loss of partons 
in the unique participant geometry of central U+U col- 
lisions pTol Il2| , and it is feasible for U+U collisions to 
extend our physics reach in a variety of observables. 

The transverse particle density increases by about 35% 
in the ideal limit of tip-to-tip U+U collisions compared 
to Au+Au at zero impact parameter. However, the prac- 
tical limitation of selecting the U+U samples of inter- 
est via conventional measures of centrality, as used in 
RHIC experiments, leads to uranium beams offering a 
smaller, yet still-worthwhile advantage. Specifically, our 
simulations suggest that the maximum achievable val- 
ues of (1/ S)dN c h/dy at RHIC could be increased from 
the present ~ 31 fm -2 to about 35 fm -2 with U+U colli- 
sions. This increase should justify the needed investment 
of effort and resources at RHIC, being of the same order 
as the measured increase for central Au+Au collisions 
when the beam energy was increased from ^/snn = 130 
GeV to 200 GeV, or just under half the increase between 
the top CERN SPS energy and y/s^N = 130 GeV. 
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